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Abstract - In 2005, Wyart et al. (Europhys. Lett., 72 (2005) 486) showed that the low frequency 
vibrational properties of jammed sphere packings can be understood in terms of a length scale, 
called £*, that diverges as the system becomes marginally unstable. Despite its tremendous success, 
it has been difficult to connect the counting argument that defines I* to other length scales that 
diverge near the jamming transition. We present an alternate derivation of I* and show that it 
arises naturally as the minimum size of a rigid cluster of particles with free boundaries. This 
understanding allows us to place a quantitative upper bound on the magnitude of f. We also 
present the first direct numerical calculation of f, and show that it is close to the upper bound. 
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O | Introduction. — Disordered solids exhibit many com- 
i^pion features, including a characteristic temperature de- 
pendence of the heat capacity and thermal conductivity [1] 
and brittle response to mechanical load [2] . A rationaliza- 
^ j t'ion for this commonality is provided by the jamming sce- 
OQnario [3,4], based on the behavior of ideal sphere packings, 
which exhibit a jamming transition with diverging length 
scales [3,5-12] as a function of packing fraction. According 
^^to the jamming scenario, these diverging length scales are 
responsible for commonality, much as a diverging length 
CO near a critical point is responsible for universality. 

I i One key diverging length scale at the jamming transi- 
!> tion arises from the so-called cutting argument introduced 
kvJby Wyart et al. [6,7]. The cutting argument is a count- 
ing argument that compares the number of constraints on 
& each particle to the number of degrees of freedom in a sys- 
tem with free boundary conditions. The cutting length, 
i* , is directly tied to the anomalous low-frequency behav- 
ior that leads to the distinctive heat capacity and thermal 
conductivity of disordered solids [1], and thus can be con- 
sidered a cornerstone of our theoretical understanding of 
the jamming transition. 

Despite its importance, however, the connection of 
the cutting length derived from counting arguments to 
physical length scales that diverge with the same expo- 
nent [5, 9, 10] has not been understood. In this paper, we 
calculate the cutting length directly and show that it has 
a simple physical interpretation as a rigidity length. This 



suggests that it might still be relevant even for systems for 
which counting arguments are less useful, such as packings 
of frictional particles [13-15] or ellipsoids [16-18]. Finally, 
we use the new interpretation of the cutting length as a 
rigidity scale to show that it is directly related to one of 
the length scales identified by Silbert et al. [5]. 

Review of the cutting argument [6, 7]. — Con- 
sider an infinite, mechanically stable packing of frictionless 
spheres in d dimensions at zero temperature and applied 
stress. Two spheres repel if they overlap, i.e. if their cen- 
ter to center distance is less than the sum of their radii, 
but do not otherwise interact. "Rattler" particles that 
have no overlaps should be removed. Since the remaining 
degrees of freedom must be constrained, the average num- 
ber of contacts on each particle, Z, must be greater than 
or equal to 2d, which is precisely the jump in the contact 
number at the jamming transition [3, 12]. All lengths will 
be given in units of <r, the average particle diameter, and 
frequencies will be given in units of y^k e g/m, where k e s 
is the average effective spring constant of all overlapping 
particles and m is the average particle mass. 

It is instructive to study a simpler system, the "un- 
stressed" system, in which each repulsive interaction be- 
tween pairs of particles in the system is replaced by a 
harmonic spring of stiffness k at its equilibrium length. 
The geometry of this spring network is identical to the 
geometry of the repulsive contacts between particles in 
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the original system and the vibrational properties of the 
two systems are closely related [7] . Now consider a square 
subsystem of linear size L obtained by removing all the 
contacts between particles (or, in the language of the un- 
stressed system, all springs) that cross the boundary be- 
tween the subsystem and the rest of the infinite system. 
Let the number of zero frequency modes in the cut sys- 
tem be q and the number of these zero modes that extend 
across the cut system be q' . Wyart et al. [6,7] used these 
modes to construct trial vibrational modes for the original 
infinite packing, as follows. If we restore the cut system 
with these q' extended zero modes back into the infinite 
system, the modes would no longer cost zero energy be- 
cause of the contacts that connect the subsystem to the 
rest of the system. Suppose we deform each mode sinu- 
soidally so that the amplitudes vanish at the boundary. 
This deformation increases the energy of each mode to 
order uj\, where ujl ~ 1/L. 

Note that if a mode is not extended, then it must be 
localized near the boundary, since the uncut system has 
no zero modes. However, the above procedure involves 
setting the mode amplitude to zero at the boundary, and 
so cannot be applied to such modes. It is therefore crucial 
to consider only the q' extended modes. 

The cutting argument now makes the assumption that 
q = aq, where a is a constant independent of L. Be- 
fore the cut, the number of extra contacts in the subsys- 
tem above the minimum required for stability is _/\Hr xtra ~ 
(Z - 2d)L d . When the cut is made, we lose 7V c cut - L d ~ x 
contacts. Naive constraint counting suggests that q' ~ q = 
max(- (iV c cxtra - 7V c cut ) ,0), as shown by the solid black 
line in fig. 1. Since 7V ( ? xtra and 2V£ ut both depend on L, 
we can define a length scale I* by 

q' = if L > I* (I) 
q'>0 iiL<£*. (2) 

The onset of zero modes is marked by J /V° xtra = N£ nt , so 



The variational argument predicts that at least q'/2 of 
the total L d eigenmodes of the full system must have fre- 
quency less than u>l , so the integral of the density of states 
from zero to ujl must be 

r L q' 
j dLoD(oj)>^. (4) 

However, D(oj) is an intrinsic property of the infinite sys- 
tem and must be independent of L. Therefore, assuming 
no additional low frequency modes beyond those predicted 
by the variational argument, we can vary L to back out 
the full density of states, as follows. 
If L > t, then q 1 = q = and 

put 

/ dujD(uj) = 0. (5) 
Jo 



For L < £*, we can write q'/2 = a(iV c cut - iV c cxtra )/2 = 
L d (wj, — l/£*), where appropriate constants have been 
absorbed into u>l and £* . This leads to 

/ dujD(uj) =uj l - l/t. (6) 
Jo 

Equations (5) and (6) imply that 

£>H = |° + l ^ <UJ l (7) 

v ' \ const. if uj > us, w 

where ui* = l/l* ~ Z — 2d defines a frequency scale. Note 
that while I* is potential independent, the units of fre- 
quency, and thus uj* , depend on potential [4]. This argu- 
ment predicts that the density of states has a plateau that 
extends down to zero frequency at the jamming transition, 
where Z — 2d = 0. Above the jamming transition, when 
Z — 2d > 0, the plateau extends down to a frequency uj* 
before vanishing. This agrees well with numerical results 
on the unstressed system [7]. Note the importance of the 
length scale £*, which defines the frequency scale uj* and 
is responsible for the excess low frequency modes. 

Need for a cleaner definition. — In the cutting ar- 
gument, the length scale £* is defined as the size of a cut 
region, L, where the number of extended zero modes, q' , 
first vanishes (eq. 2). The argument then assumes that 
this occurs when the cut system is isostatic, i.e. when 
jycut _ jyractra_ g e j ow ^ we b r i e fly ciiscuss a few complica- 
tions that arise with this definition and motivate the need 
for a cleaner, more physical understanding of £*. 

Shape of the cut. Wyart et al. only consider cutting 
along a flat surface and, while they argue that this is a 
reasonable choice for their variational argument, one can 
still imagine taking a cut of a different shape. If one were 
to consider a shape with a non-trivial fractal dimension, 
then A^ ut would no longer scale as L d_1 , resulting in a 
length scale with entirely different scaling. Without some 
physical motivation for what cut to make, the scaling of 
£* is not well defined. 

Dependence on the choice of degrees of freedom. The 
definition of £* in terms of the number of zero modes in 
eq. 2 is sensitive to the choice of degrees of freedom. For 
example, rattlers must be removed and internal degrees of 
freedom like particle rotations must be suppressed. For 
packings of ellipsoidal particles, to take one example, the 
choice of degrees of freedom is critical. Ellipsoids jam be- 
low isostaticity [16] and can have an extensive number of 
zero modes, suggesting that £* — oo because cutting con- 
tacts can only increase the number of zero modes. How- 
ever, when the aspect ratios of the ellipsoids are small, 
there is a band of modes similar to those for spheres, 
with a density of states that exhibits a plateau above 
uj* ~ Z — 2d [17]. One would expect a length scale 
£* ~ l/uj*, but constraint counting fails to predict this. 

Zero modes always exist at the boundary of a cut. The 
cutting argument assumes that the number of extended 
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zero frequency modes of the cut system, q', is proportional 
to the total number of zero modes of the cut system, q, 
with a proportionality constant a that does not depend on 
the size of the cut system, L. This assumption is critical 
to the derivation of the scaling of £*. In rcf. [7], Wyart et 
al. show numerically that this is true for isostatic systems 
where Z = 2d, but they do not provide such evidence for 
over-constrained systems. We find that for Z > 2d, zero 
modes exist for all system sizes, implying that the relation 
q' = aq does not always hold. 

To demonstrate this, we generate packings of N = 
4096 frictionless disks in 2 dimensions at zero tempera- 
ture. Particles i and j interact with a harmonic, spher- 
ically symmetric, repulsive potential given by Vfaj) = 
| (1 — rij/aij) 2 only if ry < cxy, where r%j is the centcr- 
to-center distance, cry is the sum of their radii and e = 1 
sets the energy scale. Mechanically stable packings were 
generated with periodic boundary conditions by quenching 
the total energy to a local minimum (using a combination 
of linesearch methods, Newton's method and the FIRE al- 
gorithm [19]). To fix the distance to jamming, the density 
was adjusted until a target pressure, p, was reached. Sys- 
tems were thrown out if the minimization algorithms did 
not converge. As in the cutting argument, we replace each 
packing with a geometrically equivalent unstressed spring 
network before measuring the vibrational modes. 

We create a cut system by first periodically tiling the 
square unit cell, consistent with the periodic boundary 
conditions. We then remove all particles whose center is 
outside a box of length L. By first tiling the system, we 
are able to take cuts that are larger than the unit cell, 
as well as cuts that are smaller. We have checked that 
our results are not dependent on the choice of N = 4096 
particles per unit cell. 
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Fig. 1: Number of excess zero modes as a function of the num- 
ber of excess contacts after the cut. Each data point is an 
average of configurations at constant pressure. 

Figure 1 shows q as a function of j\r? xtra — N^ nt . Evi- 
dently, q > for all cut sizes L and values of N^ xtra — N^ ut . 
From the cutting argument, this would either imply £* = 




small 
rigid 
clusters 



+ 



no rigid 
clusters 



macroscopic 

rigid 
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Fig. 2: a) An arbitrary surface (solid black line) of size L and an 
enclosed rigid cluster (stripes). The rigid cluster has an average 
contact number of Z' in the bulk and Z at the boundary. As 
L becomes large, fluctuations in Z' and Z vanish, b) Possible 
values of L such that a rigid cluster fits within the surface. 
Rigid clusters can either be small or larger than some minimum 
value. This minimum value defines £* . 



oo, or q' aq. Note that these zero modes arise from lo- 
cally under-constrained regions, despite the fact that the 
system is over constrained at the global level. 

The cutting argument can be rescued by assuming that 
the excess zero modes are associated with the boundary 
and that their number scales as L d_1 [6]. However, we find 
that these modes penetrate a non-negligible distance into 
the bulk of the system and that the scaling is not obvious. 

Cluster argument. — We now reformulate the cut- 
ting argument in a way that does not rely on the total 
number of zero modes but is specifically designed to iden- 
tify the onset of extended zero modes, which are the ones 
needed to obtain u>*. The counting will be very similar to 
that in the cutting argument, but the setup and interpre- 
tation will be different. Our argument will be motivated 
by the simple fact that if none of the zero modes are ex- 
tended, there must be a cluster of central particles that 
these modes cannot reach. Any deformation of this clus- 
ter costs energy, meaning it is rigid and has a non-zero 
bulk modulus. We will see that there is a strict lower 
bound for the size of such a rigid cluster, and that this 
minimum size defines £*. Thus, £* is the minimum size of 
a finite jammed system with free boundary conditions. 

In our reformulation, we define a rigid cluster as a group 
of particles in an infinite d dimensional system with av- 
erage contact number Z such that, if all other particles 
were removed, the only zero modes in the unstressed sys- 
tem would be those associated with global translation and 
rotation. This is purely a geometrical definition and is in- 
dependent of potential. Now, consider an arbitrary d — 1 
dimensional closed surface with characteristic size L (for 
example, the solid black curve in fig. 2a). We will now 
ask whether or not it is possible for the cluster of particles 
within this surface to be rigid. 
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For the cluster to be rigid, it must satisfy 

+ 1) 



(8) 



where N and N c are the number of particles and contacts 
in the cluster, respectively, and ^d(d + 1) is the number of 
global translations and rotations. This is a necessary but 
not sufficient condition for rigidity. We can write -/V c as 

N c = -Z'(N - N hndly ) + -ZN hndly , (9) 

where Z is the contact number of the iVbndry particles on 
the boundary and Z' is the contact number of the particles 
not on the boundary (see fig. 2a). We can also define 
the positive constants a and b such that N = 2aL d and 
-^bndry = 2bL d ~ 1+1 , where 7 > depends on the shape 
of the surface, with 7 = for non-fractal shapes. For 
shapes that have multiple characteristic lengths, e.g. a 
long rectangle, the choice of which length to identify as L 
is irrelevant as it only leads to a change in the constants 
a and b. In practice, however, we will always take L to be 
the radius of gyration. 

Equation (8) now becomes 



aL 



d-1+7 



{{Z 1 - 2d)L 



1-7 



c)>-id(d + i), (10) 



where c = -(Z' — Z) > 0. Equation (10) is trivially satis- 
fied if (Z' - 2d)L 1 -i -c>0, which implies 



(Z',c,7) = 



Z' - 2d 



1/(1-7) 



(11) 



We will refer to clusters that satisfy eq. (11) as macro- 
scopic clusters. However, it is also possible for (Z' — 
2d)L 1 ~ 1 — c < 0, provided L is very small, because the 
right hand side of eq. (10) is small and negative. 

It follows that it is only possible for the particles in 
our arbitrary surface to form a rigid cluster if the cluster 
is either very small or larger than L m i n ; rigid clusters of 
intermediate sizes cannot exist! 

Note that if L is large, then fluctuations in Z' and c 
vanish and Z' = Z. L mm is thus constant for all transla- 
tions and rotations of the surface and is independent of L, 
depending only on the actual shape of the surface. 

Given our arbitrary shape parameterized by c and 7, 
L m i n (Z, c, 7) is the minimum possible size of any macro- 
scopic rigid cluster in the Z — 2d <C 1 limit. However, 
we wish to find the minimum size of any rigid cluster re- 
gardless of shape, which we do by finding c* and 7* that 
minimize L m ; n and defining £* = L m [ n (Z, c* , 7*). In the 
limit Z — > 2d, wc immediately see that 7* = and 



(12) 



Z -2d' 



As depicted in fig. 2b, we are left with the result that rigid 
clusters must either be very small or larger than £*. 
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Fig. 4: a) I* as a function of Z — 2d, measured for individual 
systems as described in the text, b) Z « 3.25 in the limit 
Z —¥ 2d, close to the predicted bound. The solid white line in 
a) is the quantitative prediction of Eq. (13) using Z = 3.25. 
c) The projection, V{ui), of the single extended zero mode just 
below t* onto the modes of frequency uj in the uncut system, 
averaged over many realizations. 



We will now derive an upper bound for the magnitude 
of £* in the Z — > 2d limit. Since c ~ LN^ndiy/N, c is min- 
imized when the shape is a d dimensional hypcrsphere. 
We can approximate N and iVb n( j ry to be N w 4>V^ and 
iVbndry ~ 4 ) S^_ 1 , where <f> is the packing fraction and 
Vj^ and <Sj_i are the volume and surface area of a d di- 
mensional hypcrsphere with radius of gyration L. Using 
S^_ 1 /V^ = Wdd/L, where Wd is the ratio of the radius of 
gyration of a hypersphere to its radius 1 , the Z — > 2d limit 
of £* becomes 



w d d(2d - Z) 
Z- 2d 



(13) 



Equation (13) is a quantitative derivation of £* as a func- 
tion of Z that depends only on the value of Z, the average 
contact number at the boundary. 

We put an upper bound on £* by obtaining a lower 
bound for Z. First, note that any particle at the bound- 
ary of the rigid cluster cannot have d or fewer contacts. 
Removing such a particle would remove d degrees of free- 
dom and at most d constraints, and so the rigidity of the 
rest of the cluster would not be affected. Thus, Z > d + 1 
and 



< 



w d d(d - 1) 
Z -2d 



(14) 



Measuring £* from rigid clusters. — Note that not 
all surfaces larger than £* contain only a rigid cluster with- 
out non-trivial zero modes. Instead, there are almost al- 
ways zero modes that exist around the boundary in par- 
ticles that are not part of a rigid cluster. This is demon- 
strated in fig. 3, which shows progressively smaller cut 
systems in 2 dimensions. The faint blue disks are those 
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L= 61.912 




Fig. 3: Subsystems cut from a N — 4096 particle system, a) A large subsystem with q = 60 non-trivial zero modes. Only 
particles circled in blue participate in the zero modes. The solid red particles form a rigid cluster, b) A smaller subsystem with 
q — 35 zero modes, c) A subsystem obtained by removing one additional particle from the system in (b). This added a single 
additional zero mode that extends across the entire system. The largest remaining rigid cluster only contains 21 particles. The 
appearance of this extended mode and the corresponding breakup of the rigid cluster marks the length scale 
£* . d) A small system below £* with q = 33 zero modes. The largest rigid cluster contains 14 particles. 



outside of the cut subsystem, while the bright blue circled 
disks are those that participate in zero-frequency modes 
in the cut subsystem. Particles that do not participate in 
any zero mode are shown in red. These red disks form the 
largest rigid cluster that fits into the square subsystem. 

To identify the particles that do not participate in zero 
modes, we use a pebble game algorithm that was devel- 
oped in refs. [20, 21] to understand the rigidity percola- 
tion transition in bond- and site-diluted lattices. This al- 
gorithm decomposes any network into distinct rigid clus- 
ters and can also be used to calculate the number of zero 
modes. 

Figure 3 illustrates the qualitative picture presented in 
fig. 2b. In fig. 3a, there is a rigid cluster that covers about 
84% of the system. In fig. 3b, a smaller subsystem is cut 
from the same original system. In this smaller subsystem, 
the largest rigid cluster is smaller than in fig. 3a, but like- 
wise covers about 84% of the subsystem. In fig. 3c, the cut 
subsystem has been reduced by just enough so it contains 
one less particle than fig. 3b. However, the macroscopic 
rigid cluster has vanished and the largest rigid cluster in 
fig. 3c contains only 21 particles! As the size of the cut re- 
gion is decreased further (fig. 3d) , the largest rigid cluster 
remains of order 10 particles. 

We emphasize that the breakup of the rigid cluster be- 
tween figs. 3b and 3c is sharp: removing one particle com- 
pletely destroys the macroscopic rigid cluster with no in- 
termediate sized rigid clusters, as predicted above. How- 
ever, the cluster in fig. 3b is not the smallest macroscopic 
cluster because there still may be particles at the bound- 
ary of the cluster that are not essential for rigidity. To 
find the smallest macroscopic cluster, we pick a particle 
at the boundary and remove it only if doing so does not 
break apart the cluster. We do this until all particles at 
the boundary have been checked, and then measure the 
radius of gyration of the cluster, which we identify as £* . 

Figure 4a shows that £* diverges as (Z — 2d) -1 , con- 



sistent with the cutting argument and our reformulation. 
Figure 4b shows that Z w 3.25 as Z — > 2d, close to 
the theoretical lower bound of 3. The solid white line in 
fig. 4 shows the quantitative prediction from eq. (13) using 
Z = 3.25, which agrees extremely well with the data. 

According to ref. [6] , the extended zero modes of the cut 
system should be good trial modes for the low frequency 
modes of the system with periodic boundaries. Consider 
a system just below £* so that there is only one extended 
zero mode. The global translations and rotations, as well 
as the boundary zero modes, can be projected out of the 
set of zero modes by comparing them to the modes of 
the system just above £*. Figure 4c shows the projection 
of that single extended zero mode onto the dN modes of 
the full uncut system as a function of the frequency of 
the uncut modes. This mode projects most strongly onto 
the lowest frequency modes, implying that it is, in fact, a 
good trial mode from which to extract the low frequency 
behavior, as assumed [6]. 

Discussion. — We have presented an alternative defi- 
nition of the length scale £* from a cluster argument based 
on rigidity. In contrast to the cutting argument, we did 
not need to assume anything about the number of zero 
modes in a cut system and we considered all possible cuts 
simultaneously. As in the cutting argument, however, we 
did assume that spatial fluctuations in Z are negligible. 
Wyart et al. argue [7] that fluctuations in Z are negligi- 
ble in d > 2 dimensions, and that the condition of local 
force balance suppresses such fluctuations even in d = 2 
in jammed packings. We have applied our procedure from 
the previous section to bond-diluted hexagonal lattices 
where these fluctuations are not suppressed. Although 
these systems display a global rigidity transition [20, 21] 
when they have periodic boundary conditions, they do not 
exhibit an abrupt loss of rigidity at some length scale that 
could be interpreted as £* when they have free boundary 
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conditions. It remains to be seen if £* exists in this sense 
for bond-diluted 3 dimensional lattices. 

The interpretation of £* presented here still depends on 
the choice of degrees of freedom. It also relies on knowl- 
edge of the vibrational modes of the unstressed system, 
which is difficult to obtain experimentally. However, our 
result that £* marks a rigidity transition suggests that the 
clastic properties of a system could be used to measure I*. 
Such a measurement should be experimentally tractable, 
would not require knowledge of the vibrational properties, 
and would not require specification of the degrees of free- 
dom of the system. 

We note that the definition of the length scale £* as the 
size of the smallest macroscopic cluster that can be rigid, 
or equivalently, that can have nonzero elastic moduli, fol- 
lows naturally from the notion of jamming as the onset of 
rigidity, or more precisely, of nonzero elastic moduli [12]. 
The rigidity interpretation of £* makes it transparently 
clear that the cutting length £* is equivalent to the length 
scale, £l ; identified by Silbert et al. [5]. For systems with 
periodic boundaries, the anomalous modes derived from 
the zero modes swamp out sound modes at frequencies 
above oj*. Thus, the minimum wavelength of longitudinal 
sound that can be observed in the system is £l = cl/lu*, 
where cl = y/B/p is the longitudinal speed of sound, B is 
the bulk modulus, and p is the mass density of the system. 

For systems with free boundaries, rigid clusters cannot 
exist on length scales below £* , so the bulk modulus and 
speed of sound vanish. The minimum wavelength of lon- 
gitudinal sound that can be supported is therefore given 
by the minimum macroscopic cluster size, £*. From the 
scalings of B and w*, we see that £\, ~ (Z — 2d)~ 1 ~ I*. 
Our definition of £* implies that the two length scales not 
only have the same scaling but have the same physical 
meaning. 

Silbert et al. also identified a second smaller length scale 
£t from the transverse speed of sound, which depends on 
the shear modulus. For systems with free boundaries to be 
rigid, they must support both longitudinal and transverse 
sound, and so while our reasoning applies to both £^ and 
£t, £* should be the larger of the two, so that the condition 
for rigidity for a cluster of size L is L > £* = l\, . Note that 
systems with periodic boundary conditions of size L ^> £t 
are stable to infinitesimal deformations of the shape of the 
boundary [22,23]. 

Our result that rigid clusters cannot exist on length 
scales below £* appears to be consistent with results of 
Tighe [24] , as well as that of During et al. [25] , for floppy 
networks below isostaticity. There, they find that clusters 
with free boundaries replaced by pinned boundaries can- 
not be rigid for length scales above 1/\Z — 2d\. The use of 
pinning boundary particles has also been used by Mailman 
and Chakraborty [26] to calculate a point-to-set correla- 
tion length above the transition that appears to scale as 

r. 

Ideal sphere packings have the special property that the 
number of contacts in a packing with periodic boundary 



conditions is exactly isostatic at the jamming transition in 
the thermodynamic limit [3,12]. Here, we have shown that 
the number of contacts in such a system with free bound- 
ary conditions is exactly isostatic (eq. (8) is satisfied with 
a strict equality) in the cluster of size £* . This simplicity 
makes ideal sphere packings a uniquely powerful model for 
exploring the marginally jammed state. 
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